Description of superdeformed bands in light N = Z nuclei by cranked HFB method 
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Superdeformed states in light N — Z nuclei are studied by means of the self-consistent crank- 
ing calculation (i.e., the P-I-QQ model based on the cranked Hartree-Fock-Bogoliubov method). 
Analyses are given for two typical cases of superdeformed bands in the A ~ 40 mass region, that 
is, bands where backbending is absent (*''Ca) and present (^^Ar). Investigations are carried out, 
particularly for the following points: cross-shell excitations in the sd and pf shells; the role of the 
g9/2 and d5/2 orbitals; the effect of the nuclear pairing; and the interplay between triaxiality and 
band termination. 
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I. INTRODUCTION 

New "islands" of superdeformation (SD) were found, 
after nearly fifteen years of systematic search, in the nu- 
clear chart around the ~ 40 mass region (e.g., ■^^Ar [1] 
and ''^"Ca Q)- Surprisingly, these light and symmetric 
(i.e., N = Z) nuclei in the latest SD "archipelago" are 
magic and near-magic systems, whose ground states have 
spherical shape. 

These light nuclear systems with magic and near-magic 
numbers need cross-shell excitations, involving both the 
sd and pf shells, in order to produce collective degrees of 
freedom necessary for the formation of SD states. The 
corresponding shell-model space becomes inevitably very 
large. However, modern high-performance computation 
systems are quickly advancing to allow shell-model diago- 
nalization to be executed if the minimum and reasonable 
truncations are justified in the model space. 

Mean-field descriptions have a numerical advantage in 
reducing the dimension over exact diagonalization, owing 
to the ansatz for a many-body wave function (for exam- 
ple, the Slater determinant in the Hartree-Fock theory). 
The mean-field approach was exclusively applied [1] to 
the early studies of SD states known before 2000 [s*] , such 
as in the A ~ 80 (e.g., ^^Zui W), A ~ 150 (e.g., ^^^Dygg 
0), and A ~ 190 (e.g., ^'^^Rgu^ M) mass regions. The 
main reason is that nuclei in these SD archipelagoes be- 
long to heavy- and medium- weight classes. They are still 
out of reach of the shell-model diagonalization approach 
using the full model space in the relevant valence shells. 
The mean-field method has other advantages, particu- 
larly related to intuitive understanding of many-body 
systems, such as nuclear deformation and nuclear super- 
conductivity. 

At present, it is true to say that neither mean-field 
calculations nor truncated shell-model diagonalizations 
are dominantly superior to their counterpart. They are 
complementary at the moment. Many theoretical studies 
using these two approaches followed after the experimen- 
tal reports were pubhshed on '^^Ar and '^''Ca. 



Shell-model diagonalizations were performed by Cau- 
rier et al. [l|, and Poves Q, with the removal of 

the ds /2 orbital from the sd-pf shells and with the corre- 
sponding effective interaction. Recently, they attempted 
to give a consistent description of the SD and normal 
deformed states in '^''Ca [l3|- These calculations with 
the truncation were very successful in reproducing the 
experimental energy spectra. 

The first attempt through the mean-field approach was 
carried out with the cranked Nilsson model in the orig- 
inal paper [l|. In the ^°Ca paper ,2], the cranked RMF 
(relativistic mean field) method was applied. Both of the 
methods ignored the pairing correlation, so that the en- 
ergy spectra of the low-spin members (J < lOS) in the 
SD bands were not well reproduced. 

Long and Sun then applied the projected shell model 
(PSM). In this model, basis are produced through 
angular momentum projection onto the Hartree-Fock- 
Bogoliubov (HFB)_states, obtained with the P+Q-Q two- 
body interaction |ll| . A merit in this framework is that 
the pairing correlation is properly treated. As a conse- 
quence, a better agreement was obtained with the ex- 
perimental data. However, the calculation was restricted 
to an axial symmetric shape (with the deformation pa- 
rameter fixed all the way from low- to high-spin regions), 
so that the shape evolution of the system, particularly 
the triaxial degree of freedom in response to the Coriolis 
force, cannot be discussed in this model. 

Variable deformation is an important degree of free- 
dom in a rapidly rotating nucleus. For example, the 
band-termination phenomenon for the SD band is pre- 
dicted by the cranked Nilsson model, which gives a con- 
tinuous evolution in triaxiality towards the non-collective 
oblate deformation (7 = —60°) in the band limit ^. In- 
akura et al., applied the cranked Skyrme Hartree-Fock 
(HF) method, which does not restrict nuclear shape un- 
like the PSM, but the pairing correlation was ignored 
Bender, Flocard, and Hcenen analyzed the SD bands in 
the A ~ 40 region by means of the most sophisticated 
method, the generator coordinate method (GCM) with 
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the projected Skyrme HF+BCS states [Tl|. The Lipkin- 
Nogami method and particle number projection were 
apphed, so that the pairing was properly treated. Al- 
though nuclear shape can vary through the constraint on 
the quadrupole moment, only axial deformation was as- 
sumed in the calculation. In addition, the analyses were 
restricted only to low-spin states (J < 6h). This is be- 
cause the states with non-zero angular momentum were 
only kinematically created through angular momentum 
projection (without cranking). A dynamical effect orig- 
inating from shape coexistence was considered through 
the GCM, but the method underestimated the more im- 
portant dynamical effect coming from the Coriolis force, 
which plays a major role at high spin. 

The aim of this paper is thus to test another mean 
field approach, which can handle the pairing correlation, 
the Coriolis force, and the evolution of nuclear shape (in 
particular, triaxiality) in a fully self-consistent manner, 
for the full sd-pf model space. For this purpose, the SD 
bands in ^^Ar (a case with backbending) and ''"Ca (a 
case without backbending) are analyzed with the P-l-Q-Q 
model based on the HFB method in this work. 

II. THE P-HQ Q MODEL OF SELF-CONSISTENT 
CRANKING CALCULATION 

In the current framework, the Hamiltonian contains 
two terms, 

H = Hn + V. (1) 

The first term Hq represents the one-body term and 
it is the spherical Nilsson Hamiltonian. In the second 
quantization notation, it is expressed as 

-ff = ^ emalnam, (2) 

m 

where a pair of operators {al^, am) denotes fermionic op- 
erators of creation and annihilation. This part is solved 
exactly 

Ho\lpm) = GmlV'm), (3) 

and the eigenstates {|V'm) = aL|0)} (the spherical Nils- 
son states) are used as the basis in the following stages. 
Index m collectively denotes the quantum numbers in the 
Nilsson model, that is, (nljft), as well as isospin and par- 
ity. A time-reversal state of m is denoted as rh. We use 
the notation "m > 0" which means (nlj;Q, > 0). In this 
case, its time-reversal state fh corresponds to [nlj] — fi). 
The so-called Nilsson parameters for the spin-orbit and 
orbit-orbit forces (denoted as k and ^ in the standard 
notation) are taken from Refs. [l^ and [ll]. 

The second term V represents the two-body part, and 
it is the P-l-Q-Q interaction in this study, 

y^-l^xY. - E GrHPr, (4) 

fi— — 2 r— p,n 



where the first and second terms correspond to the 
particle-hole and particle-particle channels of the two- 
body interaction, respectively. The former interaction 
is responsible for the long-range correlation to describe 
nuclear deformation, while the latter is for the short- 
range correlation to handle the nuclear pairing. The 
quadrupole operator and the monopole pairing operator 
are respectively given as 

= '^iQt^)rnn"'ln^n, (5) 

Pr = ^ ama„. (6) 

m(er)>0 

The Hamiltonian is "diagonalized" with the basis ipm, 
by means of the mean field approximation. It corre- 
sponds to a procedure to extract one-body ingredients, 
Vmf, from the two-body interaction, V, so as to diagonal- 
ize Vmf- The residual part, = V — Vmf, is therefore 
neglected in the approximation. 

Remembering that we take the pairing correlation into 
account, the HFB ansatz is employed in the variational 
calculations. That is, 

|HFB)=n/3p|0). (7) 

p 

The quasiparticle annihilation-creation operators {/3p,Pjj) 
are related to the original annihilation-creation operators 
{am,o-m) through the Bogoliubov transformation, 

(,^0^(v-)(<:0- 

Matrix elements U and V in the Bogoliubov transforma- 
tion correspond to the variational parameters in the HFB 
theory. The density matrix p and the pairing tensor k 
are introduced here, as 

Pran - (HFB|at a™|HFB> - (l/*y^),„„, (9) 
Kmn = (HFB|a„a™|HFB> = (l/*C/^)„„. (10) 

The mean-field approximation of the Hamiltonian thus 
reads 

Hhfb = E-^-pB +^'^Ep(3^f3p (11) 
p 

= £'hfb + hmnalnan + ^ ^ A^„a„a„ -|- h.c. 

mn r— p,n mnfEr 

The one-body component in the particle-hole (ph) 
channel, represented as h, has the following form. 

h = e + r (12) 

The first term (e)jj = SijCi represents the spherical 
Nilsson energy. The second term F is called the self- 
consistent potential, 
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where the self-consistent coefScient is given as 

Q^ = -xTr(pQ^). (14) 

The coupUng constant x is determined in the standard 
manner by comparing to the axially deformed Nilsson 
model with deformation P'^ in the beginning of the vari- 
ational calculation. 

The one-body component in the particle-particle (pp) 
channel, denoted as A in Eq. lfTTjl . describes the pairing 
correlation. It has the form, 

K^n = \p*riPr\nn (r = P,n), (15) 

where the pairing matrix element (^V),„„ is determined 
from Eq. ^ and the self-consistent pairing coefficient Pr 
is expressed as 



m{eT)>0 



(16) 



The pairing-gap energy (A) is defined as the average of 
the matrix elements of A, that is. 



M' 

A" - — V A" - 



(17) 



m(eT)>a 



where M' = M/2 is the half the dimension of a subspace 
of given isospin (r). In the case of the present separable 
interaction, the expression for the pairing gap is simply 
given as 



(18) 



The pairing strength is determined in the standard 
manner by using the Nilsson-fBCS calculation with the 
initial values for the pairing gaps (A°) together with the 

High-spin states are produced with the self-consistent 
cranking model. That is, the variational equation, 

6{BFB\H ~ LoJ^ ~ ArA^rlHFB) = 0, (19) 



the chemical potential. The presence of the chemical po- 
tential is due to the introduction of the BCS-type pairing 
correlation, which breaks the particle number conserva- 
tion. As a result, the mean particle number needs to be 
constrained in the calculation. 

The HFB energy -Ehfb can be thus written as 



HFB 



(HFB(J)|iJ|HFB(J)) = Tr(p/i) 



r— p.n 



A 



(22) 

and this corresponds to the yrast spectrum. 

For the model space, two major shells {N = 2, 3, or the 
so-called sd-pf shell) each for protons and neutrons are 
used, that is, d5/2,Si/2,d3/2 (A^ = 2); f7/2,P3/2,f5/2,Pi/2 
{N — 3). This choice is in accordance with the Kumar- 
Baranger prescription for the P-l-Q-Q force tlTj. When a 
role of an intruder gg/2 {N — 4) orbital is discussed, it is 
also included in the model space. 

Further details of the method are available in Ref. [13] ■ 



III. NUMERICAL RESULTS 

Two N = Z nuclear systems will be studied in this pa- 
per, which are ^'^Ca and '^^Ar. The former nucleus shows 
no sign of backbending so far (up to J = 16h), while the 
latter has a clear backbending at J = lOh. With the 
P-f Q-Q model based on the cranked HFB approach, we 
attempt to describe the superdeformed states in these 
nuclei in a self-consistent manner both in the ph- and 
pp-(hh-) channels. 



A. A case of no backbending: *°Ca 

The SD band of "^''Ca is so far identified up to J = 16h 
This rotational band is regular and no backbending 
is currently observed. 



1. A role of the orbital 



is self-consistently solved by means of the gradient 
method under the following two constraints: 

(HFB I J, I HFB) = Tr (pj,) = J, (20) 

where J is the total angular momentum, and 

(HFB|Ar|HFB> = Tr (p) = N, (21) 

where N is the total particle number. 

In this study, the usual one-dimensional cranking 
model is implemented, so that only one component of the 
total angular momentum vector is constrained. Quanti- 
ties Lo and Xr in the variational equation are the Lagrange 
multipliers. The first multiplier is interpreted as the ro- 
tational frequency, while the second multiplier stands for 



As stated above, our model space contains the full 
sd- and pf-shells. It is thus possible to examine, in the 
framework of the self-consistent mean-field calculation, 
the core-polarization effect, or the influence coming from 
the d5/2 orbital truncated in the shell- model calculations. 
The occupation numbers of each single-particle orbital 
provides us useful information for this aim. 

Before analyzing our own calculations, it is worth 
learning the results obtained by others. In a description 
of the SD band of *°Ca through the shell- model diago- 
nalization by Poves Q , the single-particle model space is 
set to be Si/2,d3/2 (A^ = 2); f7/2,P3/2,f5/2,Pi/2 {N = 3). 
(This choice of the model space was also used for '^^Ar.) 
The d5/2 orbital is excluded from the valence space for 
the numerical reason. Within this model space (and the 
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corresponding effective interaction), the 8p-8h configura- 
tion, that is, (si/2d3/2)^(fp)^, was proposed for the de- 
scription of the SD band. The calculated result based 
on this configuration reproduces the experimental data 
quite well with a well-tuned effective interaction [l|, 
Long and Sun raised the question about the ds /2 trunca- 
tion in their paper where they performed the PSM (pro- 
jected shell model) calculationfll]. According to their 
analysis, the d5/2 orbital does not contribute to the SD 
state in ^^Ar, but to higher excited rotational bands. 
The cranked Nilsson calculation performed for ^^^Ar pLj 
is also informative. It was obtained that only about half 
a particle in each isospin sector is excited into higher or- 
bitals from the d5/2 orbital. Summarizing these results, 
it can be said that Poves' prescription for the truncation 
might be a good approximation for a description of the 
SD band. 

Now, let us turn to our calculation. The initial pa- 
rameters for the self-consistent iterations are chosen to 
be /30 = 0.6 and (A0,A0) = (1.360,1.513) MeV. These 
initial pairing-gap energies are employed from Ref. [iS^. 

In Fig. [1] the average occupation numbers, which are 
diagonal elements of the density matrix, p, are presented 
as a function of the total angular momentum. The right 
and left panels in the figure show the occupation numbers 
for protons and neutrons, respectively. The graphs look 
quite similar to each other as a consequence of iV = Z . 
Following the presentation in Ref. [l|, the occupation 
in the d5/2 orbital is plotted using the hole occupation 
number. Our result shows (see Table |T| that only half a 
particle is missing from the fully filled ds /2 orbital, which 
supports the shell-model truncation by Poves. This result 
is similar to the cranked Nilsson calculation for ■^^Ar 1]. 
However, the number increases at higher spin (J > IQK) 
to reach to 1 (see Table |ll]) , so that the core polarization 
may need to be taken into account, particularly at high 
spin. 



Orbital 


d5/2 


Sl/2 


d3/2 


f7/2 


P3/2 


h/2 


Pl/2 


Proton 


5.48 


1.01 


1.54 


2.63 


0.88 


0.29 


0.17 


Neutron 


5.45 


1.03 


1.56 


2.55 


0.92 


0.30 


0.19 


Total 


10.93 


2.04 


3.10 


5.18 


1.80 


0.59 


0.36 



TABLE I: Occupation numbers of '"'Ca at J = 0. The sub- 
space (si/2d3/2) is occupied by about five (= 5.14) particles, 
while the pf-shell is filled with about eight {= 7.93) parti- 
cles. The hole occupation number in the Ar-,/2 orbital is 1.07 
(= 12 - 10.93). 

Table [J displays the details of the occupation num- 
bers at J = 0. The net occupation number in the d5/2 
orbitals are 10.93, that is, about one particle (in the 
isoscalar basis) is excited into upper orbitals, as already 
mentioned above. The total numbers of the occupation in 
the subspace (si/2d3/2) and the pf shell (f7/2P3/2f5/2Pi/2) 
are about five and eight, respectively. In other words, 
our calculation suggests (d5/2)~^(d3/2Si/2)^(fp)® for the 



Orbital 


d5/2 


Sl/2 


d3/2 


f7/2 


P3/2 


f5/2 


Pl/2 


Proton 


5.11 


1.22 


1.67 


2.68 


0.57 


0.52 


0.23 


Neutron 


5.11 


1.22 


1.64 


2.64 


0.61 


0.50 


0.24 


Total 


10.22 


2.44 


3.31 


5.32 


1.18 


1.02 


0.47 


TABLE II: 


Occupation 


numbers of 


Ca at J 


= 18fi,. 


The 



subspace (si/2d3/2) is occupied by about six (= 5.75) par- 
ticles, while the pf shell is filled with about eight (= 7.99) 
particles. The hole occupation number in the d5/2 orbital is 
1.78 (= 12 - 10.22). 



band-head structure of the SD band. To a good extent, 
this configuration is consistent with the 8p-8h structure 
proposed in the shell-model calculation. Although our 
calculation indicates a possible core polarization, this ef- 
fect can be minor at low spin. 

Table |ll] shows the occupation numbers at J = 18/i. 
The number of particles in the pf shell is maintained 
to be eight (= 7.99), but about one more particle 
is excited from the d5/2 to the (si/2d3/2) subspace. 
The corresponding configuration is thus approximated as 

(d5/2)-'(d3/2Sl/2)^(fp)^ 

Comparing these two tables (and also from FiglT]), it 
can be seen that the rotational band is created through 
two modes. One is excitation within the sd shell, mainly 
an excitation from the d5/2 orbital to the upper sd shell 
(si/2d3/2); the other is within the pf shell and the rel- 
evant excitation is mainly from the P3/2 orbital to the 
f5/2 orbital. It is also learned from the calculations that 
the numbers of particles in the f7/2 and the d3/2 orbitals 
are almost constant in a wide range of the total angular 
momentum. 

From the above analysis, it can be said that the main 
part of the superdeformed structure is determined by the 
eight particles in the pf shell. Whereas, rotational mem- 
bers of the SD band are mainly produced by gradual ex- 
citations from the do /2 orbital to the upper sd shell in our 
model, in addition to the minor internal reconfiguration 
inside the pf shell. 



2. A role of the 39/2 orbital in J < 16h and rotational 
energy 

It is worth examining here an effect of the gg /2 orbital 
{N — 4), which is missing from the shell model calcula- 
tion. 

Let us see first the calculated rotational energy of the 
SD band with the original single-particle model space, 
that is, without the gg/2 orbital. The calculated and 
observed spectra of the SD band are plotted in Figl^l 
The band- head {J — 0) energy of the calculated spectrum 
is normalized with the experimental data, E(J = 0) = 
5.218 MeV. The agreement looks very good, particularly 
at low spin (J < 12h). In FiglSl the excitation energy is 
plotted, following Ref. relative to a rigid rotor energy. 
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Total Angular Momentum Total Angular Momentum 

FIG. 1: Proton (left) and neutron (right) occupation numbers as a function of tiie total angular momentum for ''"Ca. In this 
case, the gg/2 orbital is excluded from the model space. 
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FIG. 2: Calculated and observed rotational energy of the SD 
band in ^"Ca. The calculated ground-state energy (at J — 0) 
is normalized with the experimental value, E{J = 0) = 5.218 
MeV. 



Er = 0.06909 J(J + 1) [MeV] (with a line labeled as 
HFB-Np, which is performed in the absence of the gg/2 
orbital in the model space). 

Despite this good agreement, the role of the gQ/2 or- 
bital is still worth an examination because the (deformed) 
Nilsson model implies that some of the split gg/2 states 
intrude into the sd shell (TV = 2) at /3 ~ 0.6. (See, 
for example, Fig 2.21a, p. 73 in Ref.p^.) For this pur- 
pose, the g9/2 orbital is added to the model space in the 
present framework, and the calculation is repeated with 
the pairing force being unchanged. The result is plotted 
with a line denoted as HFB-Gp in Fig. [3] The low-spin 
behavior is almost identical to the previous case (HFB- 
Np), that is, without the gg/2 orbital. A small deviation 
from the HFB-Np can be seen at high spin (J > I2h), 
but in practice this difference is negligible as far as the 
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FIG. 3: Excitation energies relative to a rigid-rotor energy, 
Er = aJ{J + 1), where a = 0.06909 (MeV) for *"Ca. The 
four calculated results (HFB-Np, Nn,Gp, and Gn) are com- 
pared with the experimental data. The HFB calculations 
with (without) the gg/2 orbital are denoted as G (N). In addi- 
tion, two different sets of the initial value for the pairing are 
chosen for each case. The set p corresponds to (Ap,AjJ) = 
(1.360, 1.513) MeV, while the set n to {A°, Al) = (0.15, 0.15) 
MeV. 



rotational energy is concerned. The occupation numbers 
are also plotted in FigHl which shows only slight differ- 
ences from FigHJ It looks that the inclusion of the gg/2 
orbital gives rise to only a minor influence to the nuclear 
structure, but it turns out to be quite essential to the 
high-spin nuclear structure of ^°Ca, through the subse- 
quent analyses. We will come back to this argument in 
connection to backbending. 
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FIG. 4: Proton (left) and neutron (right) occupation numbers as a function of the total angular momentum for ^''Ca. In this 
case, the g,g/2 orbital is included in the model space. 



3. An effect of pairing correlation and backbending 

An effect of the pairing correlation can be also studied 
here. Our approach here is to compare two cases: with 
and without the pairing. The case without the pairing is 
constructed by choosing the initial pairing-gap energies 
to be smaU: (AO, A^) = (0.150,0.150) [MeV]. With this 
choice, the gap energies disappear as early as J ~ O.lh 
and remain to do so at higher spin. This calculation 
is essentially the Hartree-Fock (HF) calculation without 
the pairing, like the cranked RMF 2| and the cranked 
Skyrme HF [H]. 

The rotational energies in Figl3] denoted HFB-Nn and 
HFB-Gn correspond to the case with and without the 
§9/2 orbital, respectively (the both cases are without the 
pairing). As seen in the figure, there is no much dif- 
ference between these two cases, but the both of them 
underestimate the experimental data. When the pair- 
ing is adequately taken into account (HFB-Np and HFB- 
Gp), there is a plateau structure in the graph at low-spin 
region, which brings a better agreement to the experi- 
mental data. However, when the pairing correlations are 
absent, the plateau structure disappears and the steeper 
curves appear. The similar results were obtained in the 
other mean-field calculations neglecting the pairing cor- 
relation, such as the cranked RMF model [2] and the 
cranked Skyrme HF calculation [l^l. In the case of '^^Ar, 
the shell model calculation p| and the PSM [ll| repro- 
duce the experimental data fairly well, and the plateau 
structure is seen in these calculations. These results sug- 
gest the importance of the higher order correlations in 
the two-body interaction beyond the mean-field level. 
However, our calculation also implies that an inclusion 
of the pp-channels in the mean-field approximation, that 
is, the pairing correlations, seem to "salvage" effectively 
the important correlations that the ph-channels in the 
mean-field approximation fail to pick up. To support 



this remark, discrepancies in the excitation energy start 
to happen (which is of the order of about 1 MeV), after 
the gap energies disappear at J > 6 — 7S (See Figsl3]and 
[3 as well as the subsequent discussion in the following 
paragraph) . 
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2 4 6 8 10 12 14 16 
Total Angular Momentum 

FIG. 5: Calculated pairing-gap energies as a function of the 
total angular momentum for ^"Ca. The identification of the 
curves are the same as FiglS] 

FigUl displays the calculated pairing-gap energies for 
protons and neutrons with the original initial values for 
(AO, AO) = (1.360,1.513) MeV. In either case, the gap 
energies disappear at J > 8h, which is roughly consistent 
with the point where the deviation in the calculated exci- 
tation energy can be seen in Figl3l In finite systems, the 
pairing correlation should persist even at high spin, like 
in the Fig. 2(d) in Ref.[iH. This "coUapse" of the pairing 
is a notorious problem in the BCS-type theory applied to 
a finite system, and this "phase transition" of the pair- 
ing gap is known to be a mere artifact of the model. In 
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reality, nuclear systems should undergo a crossover, or a 
gradual decrease in the pairing gap. In the present work, 
the disappearance of the calculated gap energies is rather 
smooth and gradual, so that the qualitative feature of the 
system might be expected to be maintained, as Ring and 
Schuck explain at p. 278 in their textbook [l9|. However, 
a more elaborate treatment to keep the pairing correla- 
tion is necessary for more accurate descriptions at high 
spin, such as the Lipkin-Nogami method [20!]. 

So far, no backbending is reported in the SD band 
of ''"Ca (until J ~ 16fi-), and our calculation is con- 
sistent with this observation (with or without an in- 
clusion of the gg/2 orbital). The shell model calcula- 
tion by Poves also reproduced this result. Interestingly, 
the shell model calculation predicts the backbending at 
higher spin (J ~ 20?i) [3,]. This angular momentum corre- 
sponds to the band termination for the 8p-8h configura- 
tion, that is, (d3/2Si/2)**(f7/2)^- According to Poves, the 
nuclear structure after the backbending is constructed 
by such configurations as (d3/2Si/2)^(f7/2)^P3/2, as well 
as similar configurations allowing excitations into higher 
orbitals in the pf shells. 

In order to see clearly how backbending occurs, the so- 
called "backbending plot" is convenient. In this paper, 
the transition energy, Ej{J), is defined as 

E^{J) = E{J)-E{J-2), (23) 

where £^^(0) = 0. Alternatively, the rotational frequency 
is defined as 

uj{J) E^{J)/2h, (24) 

and w(0) = 0. 



HFB-Np 

HFB-Gp 




Rotational Frequency (MeV) 

FIG. 6: Backbending plot of the SD band of '"'Ca. HFB-Np 
(HFB-Gp) corresponds to the calculation without (with) the 
g9/2 orbital in the model space. 

In Fig[6l two cases of the calculations are plotted (with 
and without the g9/2 orbital). The line denoted by HFB- 
Np (HFB-Gp) corresponds to the case without (with) 
the g9/2 orbital in the model space. As confirmed in 



Figsd] and [H there is no much difference between these 
two cases in the spin range J < 2Qh. This situation is 
reflected in FiglS] showing completely the same behav- 
iors in the two calculated lines in J < 20?i. A difference 
happens beyond J — 20h. The case without the gg/2 
orbital (HFB-Np) shows no sign of backbending even at 
as high as J = 25h ^, while the line of HFB-Gp starts 
to backbend at J = 20h, as predicted by the shell model 
calculation. However, it should be noted that the shell 
model calculation does not contain the gg/2 and d5/2 or- 
bitals. 



4- A role of the gg/2 orbital in the backbending 

To study the rotational alignment, it is useful to cal- 
culate the single-particle angular momentum component 
along the cranking axis. The quantity is given as 

Uxim)) ^^pjnn{jx)nni, (25) 
n 

where the indices m, n denote the spherical Nilsson basis. 

Fig[7] shows the single-particle spin component along 
the cranking axis, in the case where the gg/2 orbital is 
excluded from the model space. The f7/2 orbital is the 
major contributor to the total angular momentum, and 
its contribution becomes gradually increased at higher 
spin states. About 60% of the total angular momentum 
is produced by this orbital (of protons and neutrons). 
The second main component is produced by the d5/2 or- 
bitals, and its percentage to the total angular momentum 
reaches nearly 20% at J = 25h. The contributions from 
the d3/2 and f5/2 orbitals are also significant at high spin. 

FigO shows the single-particle spin components when 
the gg/2 orbital is included in the model space. Beyond 
J = 20h, a clear structural change is seen. That is, the 
contributions from the f7/2 orbital as well as d3/2 and 
^5/2 orbitals are saturated. Instead, a rapid alignment 
of the gg/2 orbital happens. Obviously, the backbend- 
ing seen in FiglS]is caused by this structural change. As 
mentioned earlier, in the case of the shell model calcula- 
tion, backbending is initiated by the band termination, 
(d3/2Si/2)'*(f7/2)*- The higher spin states are produced 
as excitations happen from the f7/2 orbital to higher or- 
bitals in the pf shell. Comparing with our results, the 
backbending occurs from similar reasons, but not exactly 
the same. First of all, the direct factor to cause the back- 
bending is the rapid alignment in the gg/2 orbital, which 
corresponds to an excitation to higher orbitals, but this 
excitation is more drastic in the sense that three major 



^ The calculation of HFB-Np can be executed up to J ~ 25h. Be- 
yond this angular momentum, the SD structure no longer exists 
in the current framework. The cranked Skyrme HF calculation 
by Inakura, et al. also shows that the SD structure ends at 
J = 24ft. 
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FIG. 8: Proton (left) and neutron (right) single-particle spin components along the cranking axis for ""Ca. In this case, the 
g9/2 orbital is included in the model space. 



shells are involved (i.e., sd-pf-g). In addition, the satu- 
ration of the alignment in the £7/2 orbital produces the 
similar mechanism as the band termination. In fact, the 
maximum angular momentum generated by two parti- 
cles occupying the £7 /2 orbital is 6h, which is almost the 
value read in Fig[8]for the £7/2 orbital. Unlike the sim- 
ple band termination picture, the d3/2 and f5/2 orbitals 
also show the saturation, despite that the generated an- 
gular momenta are less than the maximum values. No 
more additional angular momentum is created by fur- 
ther alignments of these orbitals in the sd and pf shells 
beyond J = 2Qh (i.e., before the backbending) . In other 
words, beyond this total angular momentum, only one 
high-j orbital (g9/2) dominantly produces an additional 
angular momentum on top of the collective angular mo- 
mentum already produced by the particles in the sd and 
pf shells. This mechanism is exactly the same as the 



original backbending mechanism in the rare-earth nuclei, 
where the ii3/2 orbital (usually, of neutrons) plays the 
same role as its counterpart, that is, the gg/2 orbital. 

The calculated occupation numbers (Ta- 
bles mil and IIV|1 imply that the configuration 
changes from (d5/2)"2(d3/2Si/2)^(fp)^(g9/2)^ to 
(d5/2)^^(d3/2Si/2)^(fp)^(g9/2)^ in the backbending 
region. Essentially, this change is brought by an excita- 
tion from the d5/2 orbital to the g9/2 orbital, while the 
configurations inside the pf shell and the subspace of the 
sd shell are relatively stable throughout the whole range 
of angular momentum. This excitation from the d5/2 
orbital to the g9/2 orbital can be understood through 
the Nilsson diagram around (3 ~ 0.6. There, a low-fi 
component originating from the g9/2 orbital, that is, 
[440] 1/2, behaves like an intruder orbital coming down to 
the region near the d5/2 and d3/2 orbitals. This situation 
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Orbital d5/2 S1/2 d3/2 £7/2 P3/2 hji P1/2 g9/2 



Proton 


5.04 


1.02 


1.38 


2.59 


0.64 


0.52 


0.25 


0.56 


Neutron 


5.06 


1.02 


1.39 


2.55 


0.69 


0.50 


0.27 


0.53 


Total 


10.10 


2.04 


2.77 


5.14 


1.33 


1.02 


0.52 


1.09 



TABLE III: Occupation numbers of "'"Ca at J = 18/i in the 
case with the gg/2 orbital included in the model space. The 
subspace (si/2d3/2) is occupied by about five {= 4.81) par- 
ticles, while the pf shell is filled with about eight (= 8.01) 
particles. The hole occupation number in the d5/2 is 1.90 
(= 12 - 10.10). 



Orbital 


d5/2 


Sl/2 


d3/2 


f7/2 


P3/2 


hli 


Pl/2 


g9/2 


Proton 


4.46 


1.00 


1.39 


2.49 


0.66 


0.59 


0.26 


1.16 


Neutron 


4.65 


0.99 


1.40 


2.45 


0.71 


0.56 


0.28 


0.96 


Total 


9.11 


1.99 


2.79 


4.94 


1.37 


1.15 


0.54 


2.12 



TABLE IV; Occupation numbers of *"Ca at J = 26/i in the 
case with the gg/2 orbital included in the model space. The 
subspace (si/2d3/2) is occupied by about five (= 4.78) parti- 
cles, while the pf shell is filled with eight (= 8.00) particles. 
The hole occupation number in the d5/2 is 2.89 [= 12 — 9.11). 



implies that these three positive-parity states can jointly 
compose the nuclear many-body state when the system 
undergoes superdeformation. To produce high angular 
momentum, it is efficient to place more particles into 
the g9/2 orbital. 

It is worth studying more closely the case without the 
g9/2 orbital, in which our calculation shows no sign of 
backbending. From Tables |V] and IVIl the configuration 
at high spin looks quite stable to keep the structure of 
(d5/2)~^(d3/2Si/2)®(fp)*. The rotational members in the 
band are produced through gradual excitations from the 
d5/2 orbital to the upper sd shell as well as a minor re- 
arrangement within the pf shell. In other words, this 
structural change is adiabatic against the increment of 
the Coriolis force. This adiabatic nature of the high-spin 
nuclear structure is nothing but to mean the regularity of 
the band. In this way, no backbending is seen in this case. 
FigE] shows that the total angular momentum is created 



Orbital 


d5/2 


Sl/2 


d3/2 


hl2 


P3/2 


hl2 


Pl/2 


Proton 


5.02 


1.29 


1.70 


2.70 


0.50 


0.59 


0.21 


Neutron 


5.02 


1.29 


1.70 


2.66 


0.54 


0.56 


0.23 


Total 


10.04 


2.58 


3.40 


5.36 


1.04 


1.15 


0.44 



TABLE V: Occupation numbers of ''"Ca at J = 20ft in the 
case with the gg/2 orbital excluded from the model space. 
The subspace (si/2d3/2) is occupied by about six (= 5.98) 
particles, while the pf shell is filled with about eight (= 7.99) 
particles. The hole occupation number in the d5/2 orbital is 
1.96 (= 12 - 10.04). 



Orbital 


d5/2 


Sl/2 


d3/2 


f7/2 


P3/2 


f5/2 


Pl/2 


Proton 


4.78 


1.46 


1.77 


2.75 


0.35 


0.73 


0.17 


Neutron 


4.78 


1.46 


1.77 


2.72 


0.39 


0.71 


0.18 


Total 


9.56 


2.92 


3.54 


5.47 


0.74 


1.44 


0.35 



TABLE VI: Occupation numbers of ''"Ca at J = 26ft in the 
case with the gg/2 orbital excluded from the model space. 
The subspace (si/2d3/2) is occupied by about five (= 6.46) 
particles, while the pf shell is filled with eight (= 8.00) parti- 
cles. The hole occupation number in the d5/2 orbital is 2.44 
(= 12-9.56). 



mainly through the gradual and monotonic alignment in 
the f7/2 orbital. In addition, the contributions of the 
d3/2 and the f5/2 orbitals are seen to be non-neghgible, 
as shown in the shell model calculation. As in the previ- 
ous case, the d5/2 orbital also contributes to the creation 
of the total angular momentum. 



5. Evolution of shape 

Finally, let us examine how shape evolution happens 
in our calculation. To show quadrupole deformation, 
the Hill- Wheeler coordinates 7) [21| are used in this 
work. That is, (Qq) cx /?cos7 and (Q2) oc (3sm^/^/2. 
In the P+Q-Q model, the proportional constant car- 
rying the unit of the quadrupole moment is given as 
{hcu/hc)'^ mc^X~^ y where the harmonic oscillator energy 
is ftw = 41^1-1/3 (MeV), and ^0=197 (MeV-fm). The 
mass for a nucleon is mc^ ~ 1 (GeV) and the coupling 
constant x is given in Eq. (ll4p . 

FigE] presents the calculated triaxiality (7) and elon- 
gation (/3) for the two cases (with and without the g9/2 
orbital, which are respectively denoted as HFB-Gp and 
HFB-Np, in the figure). Until J = 20h, the elongation 
does not change significantly. The higher the total an- 
gular momentum, the more shrinkage of the deforma- 
tion can be seen along the longest principal axis of the 
quadrupole moment. But the two curves start to deviate 
from each other beyond the angular momentum J = 20h. 
The curve corresponding to the case without the gg/2 
orbital (HFB-Np) shows a monotonic decrease to reach 
/3 ~ 0.42, while the curve corresponding to the case with 
the g9/2 orbital (HFB-Gp) stops decreasing at J = 20h to 
maintain the deformation larger than /3 = 0.55. The shell 
model calculation by Poves [3] also implies the shrink of 
the shape until J = 18h in the calculation of the intrin- 
sic quadrupole moment, Qo- The cranked Skyrme-HF 
calculation by Inakura et al. [12] also shows the shrink, 
and the elongation is demonstrated to keep (3 > 0.5 (for 
J < 24fi,) for the three parameter sets (SHI, SkM*, SLy4). 
This result is consistent with our HFB-Gp case, that is, 
the case with the gg/2 orbital (and with the non- vanishing 
initial pairing-gap parameters). 

On the contrary to the elongation, there is a significant 
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FIG. 9: Calculated triaxiality 7 (left panel) and elongation (3 (right), as functions of the total angular momentum for ^"Ca. 
Captions HFB-GP and HFB-NP denote the cases with and without the gg/2 orbital, respectively. 



difference in the triaxiality between the two cases (i.e., 
with and without the gg^2 orbital). The case without the 
gg/2 orbital shows that the triaxiality is always positive 
(7 > 0°) and the triaxial deformation already starts to 
grow at low spin although the triaxiality is not so large 
(0° < 1 < 10°). The other case (with the gg/2 orbital) 
shows that the nucleus is axially symmetric until the 
backbending starts at J = 20h. Beyond the backbend- 
ing angular momentum (J — 20h), the triaxiality quickly 
starts to grow with negative values, which are consistent 
with the picture of the band termination. The amount 
of triaxiality is, however, not so substantial (I7I < 10°) 
even at J ~ 25h. The cranked Skyrme-HF calculation 
[1^ gave the consistent results with our calculation for 
the case with the gg/2 orbital, that is, the 7 is negative 
and the amount (an absolute value of 7) is less than 10° 
for J < 24:h. 



6. Summary for Ca 

Let us summarize here our analysis on the SD states 
of ^°Ca. The structure of the SD band of ""^Ca is mainly 
determined by the eight particles placed in the pf shell. 
About five to six particles sit in the (si/2d3/2) subspace of 
the sd shell, which corresponds to the configuration of six 
to seven holes. The truncation of the d5/2 orbital in the 
shell model calculation demands the eight-hole configu- 
ration in the subspace, which is approximately consistent 
with our result. However, for more accurate descriptions, 
the d5/2 and g9/2 orbitals need to be taken into account 
in the model space. This is particularly so for the de- 
scription of the nuclear structure at high spin (J > 20fi.), 
where the backbending is predicted. 

In the presence of these additional orbitals, the phe- 
nomenon similar to the band termination starts to occur 
at J = 2Qh in our model, but the saturation of align- 



ments in the pf and sd shells are assisted by the quick 
alignment in the g9/2 orbital, to which particles are ex- 
cited from the d5/2 orbital. This mechanism is consistent 
with the deformed Nilsson model, where the g9/2 orbital 
comes down to the sd shell, as an "intruder" orbital at 
superdeformation (/3 ~ 0.6). As a result of the satura- 
tion of the spin alignment in the sd-pf shell and the quick 
alignment of the g9/2 orbital, triaxial deformation starts 
occur, but its amount is not so substantial (I7I < 10°) 
that one can conclude that the nucleus keeps near-axial 
symmetry with superdeformation (/3 ~ 0.6). 



B. A case of backbending : ''^Ar 

The second targeted nucleus in this paper, ■^^Ar, shows 
a backbending at J = IQh in its superdeformed band, 
which is currently identified up to J = I6/1 pj]. The 
(si/2d3/2)'*(pf)^ structure is proposed in the shell model 
calculation truncating the ds /2 orbital [l| . For the cause 
of the backbending, simultaneous alignments of protons 
and neutrons in the f7/2 orbitals were suggested by the 

PSM [nl. 



1. Deformation of the band head 

There are slight disagreements in the calculated band- 
head deformation among different models. In the cal- 
culation with the cranked Nilsson model, minimization 
of the potential energy surface E{-^, (3) gave /3 ~ 0.45 
and 7 = 0° for the band head [l|. The PSM calcula- 
tion assumed axial symmetry (7 = 0°) and the fixed 
value f} ~ 0.48, no matter how high (or low) the total 
angular momentum is [Tlj . The calculations with the 
cranked Skyrme HF (the SIII and SkM* parameteriza- 
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tions for the interaction ^) resulted in /3 ~ 0.5 and 7 = 0° 
[l^ . In the calculation done by Bender, et al. assuming 
axial symmetry there was no SD minimum found in 
their mean-field solution (The Skyrme HF+BCS with the 
SLy6 force parameterization), but the projected solutions 
gave rise to the minimum at /3 ~ 0.5. 

Considering these calculations, the initial deformation 
parameter /3" = 0.5, as well as axial symmetry (7" = 0°), 
seems reasonable for our initial deformation parame- 
ters. The initial pairing-gap energies are employed to 
be (A°,AO) = (1.70,1.65) MeV, which are about 6% 
stronger than the values suggested in Ref.[ll|. This ad- 
justment is made so as to avoid the sudden disappearance 
of the gap energy. With this slight modification, the gap 
energy turns to disappear more gradually and smoothly 
at J ~ 7ft, as shown in FiglTOl We have confirmed, how- 
ever, that the overall qualitative nature of the SD state 
are unchanged for the modification. 



Orbital 


d5/2 


Sl/2 


d3/2 


f7/2 


P3/2 


f5/2 


Pl/2 


g9/2 


Proton 


5.07 


1.03 


1.23 


1.46 


0.61 


0.12 


0.11 


0.38 


Neutron 


5.07 


1.03 


1.23 


1.40 


0.65 


0.12 


0.12 


0.36 


Total 


10.17 


2.06 


2.46 


2.86 


1.26 


0.24 


0.23 


0.74 



TABLE VII: Occupation numbers of ^''Ar at J = in tfie 
case with the gg/2 orbital included in the model space. The 
subspace (si/2d3/2) is occupied by about four to five (= 4.52) 
particles, while the pf shell is filled with four to five (= 4.59) 
particles. The hole occupation number in the d5/2 orbital is 
1.83 (= 12-10.17). 



Orbital 


d5/2 


Sl/2 


d3/2 


f7/2 


P3/2 


f5/2 


Pl/2 


g9/2 


Proton 


5.42 


1.26 


1.17 


1.86 


0.11 


0.03 


0.00 


0.15 


Neutron 


5.43 


1.26 


1.17 


1.84 


0.12 


0.03 


0.01 


0.14 


Total 


10.85 


2.52 


2.34 


3.70 


0.23 


0.06 


0.01 


0.29 
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FIG. 10: The pairing-gap energies at the total angular mo- 
mentum < J < 16h. The g9/2 orbital is included in the 
model space. 



2. Backbending and roles of the gg /2 orbital 

To examine the quality of our calculation, it is useful 
to see the backbending plot first, which is presented in 
FigEH Due to the disappearance of the pairing gap at 
J ~ 7h, backbending starts earlier (at J ~ 7h) than the 
experiment (at J ~ lOh). Having accepted this discrep- 
ancy, the calculation manages to reproduce the qualita- 
tive behavior of the backbending profile of the SD band 
in this nucleus. Hence, in order to discuss the structural 
change causing the backbending, it is sufficient to study 



^ The SLy4 parameter results in triaxial deformation for the SD 
bandhead, but this result may be "less reliable" than the other 
parameter sets according to the authors Il2| . 



TABLE VIII: Occupation numbers of ^''Ar at J = 16ft in 
the case with the g9/2 orbital included in the model space. 
The subspace (si/2d3/2) is occupied by about five {— 4.86) 
particles, while the pf shell is filled with four {— 4.00) parti- 
cles. The hole occupation number in the d5/2 orbital is 1.13 
(= 12 - 10.85). 

the configurations before and after J ~ 7ft in our model. 
Let us see the corresponding occupation numbers next. 




Rotational Frequency (MeV) 

FIG. 11: Backbending plot of the SD band of ^''Ar. HFB-Np 
(HFB-Gp) corresponds to the calculation without (with) the 
g9/2 orbital in the model space. 

The occupation numbers are displayed in Fig ll2[ as 
well as in Tables IVlIl and IVlIIl Table ETH which presents 
the configurations before the backbending, shows that 
approximately four and a half particles occupy the sub- 
space of the sd shell, as well as the pf shell. The gg/2 
orbital is occupied by less than one particle, so that the 
orbital is expected to play only a limited role in the low- 
spin structure of '^^Ar. One can notice that there are 
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Total Angular Momentum Total Angular Momentum 

FIG. 12: Proton (left) and neutron (right) occupation numbers as a function of the total angular momentum for '^^Ar. In this 
case, the gg/2 orbital is included in the model space. 



about two holes in the orbital. This result im- 

plies the structure of (d5/2)"^(si/2d3/2)^-^(fp)'^-^(g9/2)^ 
at the band head of the SD band, which is similar to the 
(si/2d3/2)^(fp)* structure suggested by the shell model 
calculation. This structure, which has a configuration 
of two holes in the d5/2 orbital and one particle in the 
gg /2 orbital, is also similar to our result for the bandhead 
structure of the SD band of *°Ca with the gg/2 orbital 
included in the model space. 

Turning a focus onto the higher-spin states inside the 
SD band, one can tell that the configuration actually 
turns out to be closer to the configuration suggested 
by the shell-model calculation. From the Table IVIIH 
our calculation suggests the (d5/2)~^(si/2d3/2)^(fp)'* at 
J — IQh. The g9/2 orbital seems to play, again, no role 
in this case, so that the state resembles the 4p-4h con- 
figuration (si/2d3/2)^(fp)* in the truncated shell-model 
diagonalization. From this analysis, we can say that the 
basic configuration does not change much between before 
and after the backbending. 

FigdHalso suggests that the gg/2 orbital is not a key 
player for the backbending. Its occupation number is 
quite low and fairly regular throughout the whole range 
of the total angular momentum. Instead, the occupation 
number in the f7/2 orbital suddenly starts to increase 
in the post-backbending region (J > 7fi) to indicate a 
rotational alignment in this orbital. This result is con- 
sistent with the interpretation of the backbending by the 
PSM [11]. The Si/2 orbital also shows the increase, but 
it is more gradual in comparison to the f7/2 orbital. The 
sharp drop of the hole occupation number in the d5/2 
orbital can be noticed clearly in the figure. This means 
that the d5/2 orbital starts to be filled by the deexci- 
tations from the upper shells (about a half particle from 
the gg/2 orbital and another half from the pf shells). This 
result implies that the deformation becomes less substan- 
tial and the cross-shell excitation is suppressed to some 



extent. We will take a closer look at a relation between 
the deformation (/3) and the position of the gg/2 orbital 
intruding into the sd shell, soon below. 

An indication of the decrease in elongation at high 
spin was reported by other calculations such as the shell- 
model diagonalization the cranked Skyrme HF cal- 
culation [l3|, and the cranked Nilsson calculation [2^ . 
Fig ll3l shows the calculated quadrupole deformation in 
terms of triaxiality (7) and elongation (/3), with our 
model. With or without the gg/2 orbital in the model 
space, there is no much difference in the manner of the 
shape evolution, as already mentioned above. Before the 
backbending, the system maintains an axially symmet- 
ric shape, while triaxial deformation starts to grow after 
the backbending (in particular, J > lOfi.) although the 
absolute value of 7 is less than 10°. The small I7I value 
explains the reason for the successful descriptions of the 
PSM calculation [ll|, as well as the projected GCM cal- 
culation p^ . 

Based on the above discussions, it seems that we can 
conclude that the g9/2 orbital plays no significant role, 
unlike the "^''Ca case. This can be qualitatively explained 
from two things. One is the position of the Fermi levels, 
which are lower than the one of the ^''Ca. That is, cross- 
shell excitations cost energetically more in "^^Ar than in 
^"Ca. The other is smaller (3 value of the quadrupole 
deformation. From the Nilsson diagram, it can be learned 
that the gg/2 orbital can behave like an intruder orbital 
when the deformation is as large as (3 ~ 0.6. However, 
when the deformation becomes smaller such as /3 ~ 0.4 — 
0.5, the gg/2 orbital remains to be high in energy and the 
mixture of this orbital to the other orbitals in the sd shell 
becomes less likely. 

If the backbending plot, Fig[Tll is carefully studied, 
however, one can find that there is a little discrepancy in 
the high-spin behavior of the curves. The second back- 
bending seems to happen at J ~ \Ah in the case without 
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FIG. 13: Calculated triaxiality 7 (left panel) and elongation /3 



Orbital 


d5/2 


Sl/2 


d3/2 


h/2 


P3/2 


k/2 


Pl/2 


Proton 
Neutron 


5.38 
5.38 


1.33 
1.33 


1.29 
1.29 


1.95 
1.95 


0.00 
0.00 


0.05 
0.05 


0.00 
0.00 


Total 


10.76 


2.66 


2.58 


3.90 


0.00 


0.10 


0.00 


TABLE IX 


Occupation 


numbers of 


Ar at J 


= 16h 


in the 



case with the gg/2 orbital excluded from the model space. 
The subspace (si/2d3/2) is occupied by about five (= 5.24) 
particles, while the pf-shell is filled with four {— 4.00) parti- 
cles. The hole occupation number in the dc,p excluded is 1.24 
(= 12 - 10.76). 



the gg/2 orbital (HFB-Np). Obviously, this result is in- 
consistent with the experimental data, showing no sign of 
the second backbending. The case with the gg/2 orbital 
(HFB-Gp) shows the closer results to the experiment. 
Apparently from the analyses of the occupation numbers 
and the deformation evolution, the gg/2 orbital plays no 
active role in the SD structure. But it seems that these 
results does not mean that one can exclude the gg/2 or- 
bital from the model space. 

Let us examine how this discrepancy happens. First of 
all, let us compare the occupation numbers in these two 
cases (i.e., HFB-Gp and HFB-Np) at J = 16h. The case 
with the gg/2 orbital (HFB-Gp) is already shown in Table 
IVIIIl whereas the case without the gg /2 orbital (HFB-Np) 
is presented in Table HXl No significant difference can be 
seen from these two Tables. Only a tiny difference can 
be seen that the concentration onto the f7/2 orbital is 
slightly higher in the latter case than the former. 

It is then worth a look at the single-particle spin com- 
ponent along the cranking axis. There are differences 
between the two cases, which are found in the f7/2 and 
d5/2 orbitals at J > 14ft. Their behaviors are plotted in 
Figim Obviously, the gg/2 orbital is not involved in the 
production of angular momentum. The low occupation 
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FIG. 14: Single-particle spin components along the cranking 
axis. Only the 17/2 and d5/2 orbitals for protons are plotted 
to compare the two cases (with and without the gg/2 in the 
model space). 



number in this orbital implies the consistency with the 
case without the gg/2 orbital. On the other hand, f7/2 
and d5/2 orbitals show slightly different behaviors in the 
two cases (HFB-Gp and -Np in the figure, which corre- 
spond to the case with and without the gg /2 orbital in the 
model space). When the gg/2 orbital is removed from the 
model space, the alignment in the f7/2 orbital is slightly 
accelerated beyond J = \Ah. Whereas, in the other case, 
the alignment of the f7/2 orbital slows down a little. The 
opposite behavior is seen in the d5/2 orbitals. From this 
analysis, it can be said that the further alignment in the 
f7 /2 orbital causes the second backbending when the gg /2 
orbital is absent. This backbending should be regarded 
as an artifact and it is caused by the secondary effect 
due to the lack of the gg /2 orbital, which creates subtle 
differences from the case with the gg/2 orbital. 
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FIG. 15: Calculated and observed rotational energies of the 
SD band in ^®Ar. The gg/2 orbital is taken into account in the 
model space. The calculated ground-state energy (at J = 0) 
is normalized with the experimental value, £( J = 0) = 4.3291 
MeV. 

Therefore, despite the discrepancy in Fig llll our pre- 
vious conclusion, that is, the gg/2 orbital plays no major 
role in the SD band of '^^Ar, still holds because the qual- 
itative characters of the SD structure does not change 
whether the g9/2 orbital is considered or not. It is true, 
however, that a better description demands an inclusion 
of the gg/2 orbital in the model space, particularly for 
the backbending plot. 
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FIG. 16: Excitation energies relative to a rigid-rotor energy, 
Er = a J {J -f 1), where a = 0.08234 (MeV), for ^«Ar. 



3. Pairing correlation and band termination 

Let us now compare the energy spectrum between the 
calculation and the experimental data, which are shown 
in Figslin] and [111 Discrepancies are seen in the back- 



bending region {8h < J < 12h). As we have discussed 
earlier, these discrepancies are the result of the pairing 
collapse (See also FiglTUl). But the qualitative behavior 
is managed to be reproduced. In particular, the plateau 
structure at low spin and a sharp slope at high spin 
are well reproduced in Fig[Tni In comparison with the 
cranked Nilsson calculation in Ref.[l|, our result shows 
an improvement as a mean-field approach, which takes 
into account the pairing. 

Without the pairing, the cranked Nilsson calculation 
suggests an occurrence of the band-termination phe- 
nomenon for the SD band of Ar [22[ . The typical fea- 
ture of the band termination appears as a shape change 
into oblate deformation. Earlier in the present study, we 
have learned that the triaxial deformation is suppressed 
with a proper treatment of the pairing (Fig |13p . It is 
interesting to study if the band termination and the as- 
sociated oblate deformation can emerge in the absence of 
the pairing in our model. 

By choosing the initial pairing-gap parameters to be 
very weak, that is, (A|5,A°) = (0.15,0.15) MeV, the 
effect of the pairing correlation is examined in the fol- 
lowing calculation. (The g9/2 orbital is included in this 
analysis.) This choice of the initial parameters causes 
the breakdown of the pairing as early as at J ~ 1.5h. Ef- 
fectively, the calculation turns to be the HF calculation 
beyond this total angular momentum, which can be com- 
parable with the cranked Nilsson calculation without the 
pairing. The corresponding cranking calculation gives 
regular solutions until J = 16ft, the band termination 
point. But beyond the termination point, the solution 
becomes irregular and shows unphysical behavior, so that 
we ignore the calculations beyond the band termination 
point in this analysis. 

In the band termination spin (J — 16ft), the occupa- 
tion numbers of the f7/2 orbital are calculated to be 1.94 
each for neutrons and protons. Therefore, the net par- 
ticles occupying the f7/2 orbital are four (i.e., = 2-1-2), 
and the result is consistent with the shell model con- 
figuration. Two particles occupying the f7/2 orbital can 
generate the maximum angular momentum of J = 6ft, 
which corresponds to the band termination. 

In FigfTTl the calculated single-particle alignments are 
plotted for protons and neutrons. Our calculation (with- 
out the paring) reproduces the above situation: the 
single-particle angular momentum carried by the ij /2 or- 
bital is almost 12ft, consisting of 6ft for protons and 6ft 
for neutrons. 

In FigllHl the calculated quadrupole deformation, that 
is, triaxiality (7) and elongation (/3), are plotted. The 
profile for the elongation does not change from the pre- 
vious case with the pairing correlation. However, the be- 
havior of the triaxial evolution is different, especially be- 
yond J = 12ft. At the band termination point, the value 
of 7 reaches —60°, indicating the non-collective oblate 
shape. From this analysis relying on the self-consistent 
cranking calculation, it is confirmed that the band termi- 
nation phenomenon happens when the pairing correlation 



15 







• *7/2 


36 a 










.• "^3/2 






.. ^5/2 

- '5/2'Sl/2 



5 10 
Total Angular Momentum 



15 P3/2'Pl/2 



Q. 

I/) 
0) 

u 
C 

IS 

a. 
i 









36 « 










.•• "^3/2 






... J^sa 

- '5/2'Sl/2 



5 10 
Total Angular Momentum 



15 P3/2'Pl/2 



FIG. 17: Proton (left) and neutron (right) single-particle spin components along the cranking axis for ^**Ar. The gg/2 orbital 
is included in the model space. 
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FIG. 18: Calculated triaxiality 7 (left panel) and elongation (3 (right), as functions of the total angular momentum for ^®Ar, 
when the pairing correlation is switched off. 



is absent from the system (or very weak). 



4- Summary for ''^ Ar 



degree of triaxiality is not so substantial that the SD 
states are well described as an axially symmetric nuclear 
many-body system. 



As a partial summary for ^^Ar, we can conclude 
that the previously proposed structure, (si/2d3/2)^(fp)'', 
seems to be a good approximation for the structure of 
■^^Ar, according to our results. It is suggested that back- 
bending in the SD band of ■^^Ar is caused by the align- 
ment in the 17/2 orbitals, as concluded by the PSM calcu- 
lation [TT| . The gg /2 orbital does not play any significant 
role. From the successful description of the SD band 
at low spin and the failure in the backbending region, 
it was demonstrated that the pairing correlation is very 
important to describe the structure of the SD band. Tri- 
axial deformation starts to occur at high spin, but the 



IV. CONCLUSIONS 

The self-consistent cranking calculation based on the 
HFB method was applied to the superdeformed bands 
of two N — Z nuclei, ""^Ca and ^^^Ar. Our microscopic 
calculations with the P-f Q-Q interaction can manage to 
give good qualitative explanations (occasionally quanti- 
tatively) to the energy spectrum, rotational alignment, 
and backbending phenomenon of these nuclear systems. 

Special attentions were paid to the roles of (1) the d5/2 
orbital, which was removed from the sd-pf model space 
in the shell-model diagonalizations; and (2) the g9/2 or- 
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bital, which belongs to a higher shell (TV = 4) than the 
sd-pf shell. The effect of the pairing correlation was also 
investigated in connection to the evolution of triaxial de- 
formation and the band termination phenomenon. 

Inside the framework of our model, it was found that 
the truncation of the d5/2 orbital can be justified as far 
as lower-spin states are considered. Whereas, high-spin 
states are found to be produced due to a gradual excita- 
tion from the d5/2 orbital to the upper sd shell. However, 
in either case of ^°Ca and ^^Ar, an inclusion of the d5/2 
orbital does not affect the nuclear structure of the SD 
states very much. 

On the contrary, the gg/2 orbital was found to change 
the nuclear structure drastically for ^°Ca: backbending 
may happen at ,/ ~ 20h. However, the orbital plays no 
significant role for '^^Ar. These differences come from 
the location of the Fermi levels and the deformation 
which is ~ 0.6 for ^^Ca while ~ 0.4 for ^^Av. This dif- 
ference infiuences the position of the g9/2 orbital as an 
"intruder orbital" into the sd shell, in terms of the de- 
formed Nilsson model. 

The pairing correlation was found to be important to 
produce a proper energy spectrum and tend to act as a 
suppressor of triaxial deformation. Without the pair- 
ing, triaxial deformation would be enhanced and the 
non-collective oblate shape (7 = —60°) would ultimately 



emerge at high spin. However, in our model taking the 
pairing correlation into account, it is observed that tri- 
axiality is suppressed to I7I < 10°. In this sense, the 
SD states of both '^''Ca and ^^Ar are nearly axial sym- 
metric in our model, which can justify other calculations 
assuming the axial symmetry. 

Despite the success for the qualitative explanations, a 
problem was recognized in relation to the pairing collapse 
at high spin. An improvement is surely necessary for 
more accurate and quantitative descriptions of the high- 
spin structure, especially around the backbending region. 

Nevertheless, through this work, the P-l-Q-Q model 
based on the cranked HFB approach was demonstrated 
to be a practical and effective model to describe high-spin 
nuclear structure showing superdeformation. 
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